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Abstract 

The extremal t process was proposed in the literature for modelling spatial extremes 
within a copula framework based on the extreme value limit of elliptical t distribu- 
tions (Davison, Padoan and Ribatet (2012)). A main drawback of this max-stablc 
model was the lack of a spectral construction such that direct simulation was infea- 
r-pl ■ sible. In this work, we propose a spectral construction for the extremal t process 

that renders direct simulation possible. It further allows to identify the extremal 
Gaussian process introduced by Schlather (2002) as a special case of the extremal 
t process. All results naturally also hold within the multivariate domain. 
> • 
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^sO . 1. Introduction 



(N 1 iDavison et all (2012) survey the statistical modelling of spatial extremes and 

provide a global view on available models and their interconnections. Among these 
t 1 models, the extremal t process represents a max-stable process that generalizes 

the t extreme value copula to infinite dimension. It is well defined, yet no direct 
construction was known back then which lead the authors to class it among copula 
models characterized by their motivation from multivariate considerations. In that 
paper, it was further illustrated in what way the extre mal t process can b e con sidered 



as a generalization of the Brown-Resnick process (|Kabluchko et all l|2009| )), and 
an application to Swiss rainfall data bears witness of its versatility for extremal 
dependence modelling. In the following, we will see that the extremal t process 
provides a natural connection between two pr ominent max-sta ble model classes, 
namely Schlather's extremal Gaussian process ( Schlather ( 20021 )) and the Brown- 



Resnick process. We conceive a spectral construction that generalizes the one of 
the extremal Gaussian process. It renders direct simulation possible for moderately 
large general degrees of freedom. 

In the remainder of this introduction, we carry together some elementary defini- 
tions and results. Section 2 establishes and recalls useful results on convergence and 
extremal dependence structures, before we present the new spectral construction in 
Section 3 and conclude with a discussion in Section 4. 

1.1. Notational conventions 

If not stated otherwise, operations on vectorial arguments like maxima or arith- 
metic operations must be interpreted componentwise in the following presentation. 
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Vectors are typeset in bold face, in particular the vector constants = (0, ...,0) T 
and 1 = (1,...,1) T . Rectangular bounded or unbounded sets are given according 
to notations like [u, v] = [m, vi] x • • • x [u<j,Wrf] or (0, oo) = (0, oo) x ••• x (0, oo). 
The complementary set of a set B in M d is written B c . The truncation operator 
x + = max(a:,0) maps negative values to 0. The indicator function of a set B is 
denoted by Xb(-)- 



1.2. Some general extreme value theory 

For a more detailled account of ma x- stability and extrem e v alue theory in general 
we ref er the reader to the textbooks of lBeirlant et al.l ([20041 ) and lde Haan and Ferreira 
(l2006h . 



1.2.1. Max-stability 

Let Z, Zi, Z2, ... be a sequence of independent and identically distributed (iid) 
random vectors in M (d > 1) with nondegenerate univariate marginal distributions. 
We say that Z follows a max-stable distribution if sequences of normalizing vectors 
a n > and b n (n = 1, 2, ...) exist such that the equality in distribution 

max an-^Zi - b n ) = Z (1) 

i— l,...,n 

holds for the componentwise maximum. A full characterization of multivariate max- 
stable distributions leads to rather technical expressions. For our purposes, it is con- 
venient to focus on the case where the distribution G of Z has common a-Frechet 
marginal distributions Gj(zj) = & a (zj) = exp(— zJ a )\(o .00) ( z j) {j = 1, d) for 
some tail index a > 0. Monotone and parametric marginal transformations al- 
low to reconstruct all admissible univariate max-stable marginal scales in ((T|) from 
this particular marginal scale. More precisely, the class of univariate max-stable 
distributions is partitioned into the class of a-Frechet distribution under strictly 
increasing linear transformations and further the so-called Gumbel and Weibull 
classes. 

With a-Frechet marginal distributions, the standard exponent measure M is de- 
fined on [0,oc) \{0} by M((0,z] c ) = -logP(Z" < z) with the convention - log = 
00 and characterizes the dependence structure in G on a standardized scale; it is 
uniquely defined by the dependence function M(z) = M((0, z] c ) which takes the 
value 00 whenever min^ Zj — such that z (jL (0, 00). Valid dependence functions 
are characterized by Af(ej) = 1 for j = 1, d, ej = (00, 00, 1, 00, oo) T a vec- 
tor with j-th component 1, and by the homogeneity property M(rz) = r~ 1 M(z) for 
r > 0. In the independence case when G is the product of its univariate marginal 
distributions, we observe M(z) = J2j=i z J l ■ The extremal coefficient M(l) € [1,(2] 
can serve as an indicator of the strength of extremal dependence, ranging from full 
dependence associated with the value 1 to independence associated with the value 
d. 

In the infinite-dimensional domain, we call a stochastic process Z = {Z(s),s 6 
S C MP} (p > 1) with a non-empty Borel set S max-stable if its finite-dimensional 
distributions are max-stable. If Zi,Z2,... are iid copies of Z, then sequences of 
functions a n (s) > and b n (s) (n > 1) exist such that {maxt^ . ^ a n (s)~ 1 (Zi{s) — 

b n (s))} £ {Z(s)}. 



1.2.2. Domain of attraction 

Let X, Xi,X2, ... be a sequence of iid random vectors in M d with distribution 
function F. For suitably chosen normalizing sequences, relation ((T|) can hold asymp- 
totically in the sense of distributional convergence with nondegenerate marginal 
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distributions in the limit Z: 

max a n _1 (Xi — b n ) -4- Z (n — > oo) . (2) 

i— l,...,n 

Then we say that the distribution F of X is in the max-domain of attraction (MDA) 
of the max-stable distribution G of Z; or we simply say that X is in the MDA of 
Z. Normalizing sequences are not unique and the limit distribution G is unique 
up to a linear transformation. If normalizing constants can be chosen such that 
all the univariate marginal distributions Gj are of the same a-Frechet type, then 
the particular choice of b n = is admissible. In this case, the convergence in 
distribution ([2]) is equivalent to 

nP(a n - 1 X £ z) ->■ M(z a ) for all z e (0, oo) . (3) 

For d = 1, we have nPfo^ 1 ! > 2) — > z~ a (z > 0), and then X is said to be 
regularly varying at 00 with index a > or just regularly varying in the remainder 
of this paper, denoted as I £ RV Q . The normalizing sequence can be chosen as 
a n = inf{ir : P(X > x) < n" 1 }. 

For stochastic processes, the notion of MDA is defined in the sense of the con- 
vergence of all finite-dimensional distributions according to © . 



1.2.3. A spectral construction for max-stable processes 

The commonly used models for m ax-stable process es are generated with the 
help of so-called spectral constructions. ISchlatherl (|2002h proposes to use a Poisson 
process {Vi} ~ PRM(v~ 2 dv) on (0, 00) and iid replicates Qi of an integrable random 
process Q, independent of {Vi} and with EQ + (s) = 1 (s € S), in order to construct 
the max-stable process 



Z = {Z(s)} = 



max ViQi(s) 

i=l,2,... 



(s e S) 



(4) 



with univariate marginal distributions of type $1. In particular, as a consequence 
of the infinite number of points in the Poisson process {Vi}, it is possible to replace 
Qi(s) by the zero-truncated value Q^(s) in this construction. Subsequently, with- 
out loss of generality we assume that the points Vi are in descending order such that 
Vx > V2 > ••■ and V\ ~ $1. We obtain extremal Gaussian processes by choosing a 
centered and appro priately scaled Gauss ian process W for Q. Other choices of Q 
were considered (cf . iDavison et aL ( 2012)), lea d ing for instance to so-ca lled Brown- 
Resnick processes (|Brown and Resnickl (|l977h ; iKabluchko et al.l (|2009t )). The de- 
pendence function of Z for a finite number of points s%, G S is 



M Sl! ... iS(i (z) = E max 



' x Q + (»i)) 



(5) 



L5. Elliptical distributions 

Definition 1.1 (Elliptically distributed random vectors). ^4 random vector 
X m R d is said to follow a (non-singular) elliptical distribution if it allows for 

a stochastic representation X = fi + R^AXJ with a deterministic location vector 
fi, an invertible d x d matrix A that defines the dispersion matrix S = AA T = 
( a jij2)i<ji,j2<d an d 0- nondegenerate random variable Rd > independent from a 
random vector U uniformly distributed on the Euclidean unit sphere Sd-i = {z G 
M. d I z T z = 1}. We call Rd the radial variable. 

Contrary to more general definitions (cf. lAnderson and Fane! (1990)), we prefer to 
remain within the framework of quadratic and non-singular A to avoid an overly 
technical presentation for the more general cases. A prominent example are elliptical 
multivariate t distributions. 
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Definition 1.2 (Multivariate t distribution). We say that an elliptically dis- 
tributed random vector X in M. d follows the multivariate t distribution with v > 
(general) degrees of freedom if d _1 R^ ~ Fd,v, where Fd >v is the F- distribution with 
density function 

x i y [xBeta(0.5d,0.5i/)]" 1 y / (dx + iy)- ( - d+l/ \dx) d ^ ■ X[o,oo)(?) 

and Beta refers to the Beta- function. We write X ~ t v (fj,, E) and P(X < x) = 
i„(x | 

As tends to infinity, the limit of the multivariate t distribution is the multivari- 
ate normal distribution with R 2 d ~ \\- The multivariate t distribution can be 
constructed as a variance mixture of the multivariate normal distribution: With 
vY^ 1 ~ Gamma(0.5z^, 2) (y > 0) and a multivariate normal random vector W ~ 
iV(0, S) that is independent of Y , we obtain y + v^W ~ t v {p-i E), cf. for in- 
stance iDemarta and McNeill (|2005h . This construction is easily generalized to the 



infinite-dimensional setting on a domain S: If W is a centered Gaussian process 
with domain S and covariance function Cov and vY^ 1 ~ Gamma(0.5i/, 2) indepen- 
dent of W, then we call the random process V^W a (centered) t random process 
on S wh ich is characterized by the degree of freedom v and the dispersion function 
Cov (cf. iRoislien and Omrd (|2006h ). 



1.4- The t max-domain of attraction 

The multivariate t distribution fulfills the MDA condition ©. For normal- 
izing constants b n = and a n = n l / v (a^ 2 cV v )j=i, j, with c v = T(0.5(p + 
l))- 1 i/-°^+ 1 0Fr(O.5^) (cf. Table 2.1 on page 59 in iBeirlant etaD (|2004T )). we 



obtain i/-Frechet marginal distributions in the max-stable limit distribution G of 
th e multivariate t distribution t v (y, E). The dependence function of G was derived 
bv lNikoloulopoulos et al.l (|2Q09f ) : Denote by E* = (cr? j 2 )j 1 j 2 tne correlation matrix 
that corresponds to the dispersion matrix E, and by • _j = (cr^ j 2 )j 1 ^jj 2 ^j , = 
E-jj — j 2 )j x =j,j 2 jtj submatrices obtained by removing some of the rows or 
columns. Similarly for vectors, we write z^j = (z%, Zj—\, Zj+i, Zd) T ■ Then 

M„, s .(z) = "£zrH v+1 ((m-i/zj)"- 1 | E*^, + l)- 1 (E^j - ■ 
i=i 

(6) 

We refer to the max-stable limit as extremal t distribution, and we call extremal 
t process the max-stable limit of a £ random process which is a generalization of 
the extremal t distribution. Its dependence structure is characterized by the set 
of dependence functions for all finite-dimensional distributions, hence by a general 
degree of freedom v > and the correlation function 

Cov* (sj, , s h ) = yj Cov( Sil , Sjl )Cov(s h , s h ) Cav(s h , s n ) {s n , s h G S) 

which corresponds to the dispersion function Cov and determines the matrices E* 
for all finite-dimensional distributions. Similar to the copula approach in the mul- 
tivariate domain, it is convenient to call extremal t process any max-stable process 
whose set of dependence functions is the same. Put differently, we allow for uni- 
variate max-stable marginal scales beyond the ^-Frechet scale. 

2. A Breiman-type result and its ramifications 



Breiman's theorem ( Breiman (Il965h : ICline and Samorodnitskvl (jl9941) ) charac 



terizes the tail behavior of the product of a non-negative regularly varying random 
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variable with a "lighter-tailed" non-negative random variable. We state a general- 
ization adapted to our purposes. 

Proposition 2.1. Let R € RV Q (a > 0) be a non-negative random variable such 
that nV(a^ n 1 R > r) ~ r~ a for r > as n — > oo. Suppose X > is a random vector 
inM. d , independent of R, such £/ia£E(X Q ) = 1. Assume further that E(X a+s ) < oo 
for some S > or that ¥(R > r) ~ cr~ Q /or some c > 0. TTien i/ie distribution of the 
vector Y = i?X fulfills the MDA condition (|2|) wi£/i a max-stable limit distribution 
G. The normalizing vector sequences a n = (a„, a n ) T and bn = in ([2]) Zearf to 
marginal distributions $ Q in G. XTie dependence function of G is 



M Y (z) 



E 



max frt^Xf) 



(ze (O.oo)) 



Proof. We observe nP(a n X PX ^ z) = nPf^max^i^..^^ 1 X J -) > a„) for z e 
(0, oo). We can apply the clas sical Breiman theorem if E[X" +(5 1 < oo and otherwise 
a variant from Lemma 2.3 in Davis and Mikoschl ( 20081) if P(i? > r) ~ cr~ a , and 
we obtain 



nP(i? max (z^Xj) > a n ) - nP(i? > a„)E 

as n -} oo. Since nP(i? > a n ) -4 1 as n - 
marginal distributions and My(z) = E[max J= i 



max (z - 1 Xj 



(ze (o,oo)) 



oo, we can deduce the a-Frechet 



We notice that the dependence function depends on R only through the index of 
regular varia tion a. The following corollary can also be deduced from Lemma 3.1 in 
Segers ( 2012 ) or from some related derivations in Remark 3.6 in lDavis and Mikosch 
(120081) : for completeness of the presentation, here we provide a proof based on 
Proposition 2.1. 

Corollary 2.1. The process {ViQi(s)} (s G S) in the spectral construction (U) 
is in the MDA of the process Z in with normalizing functions a n (s) — n and 
b n (s) = (n > I). 

Proof. It is easily seen that the random vector V^i(Qi(si), Qi{sd)) T and its zero- 
truncated version Y = V\ {Qf(si ),..., Qi(sd)) T belong to the same MDA (if any) 
for any finite collection of points sj (j = 1, d). Hence the statement follows from 



a straightforward application of Proposition 12.11 with a = 1, R = V% ~ $i and 
X = (Qi(s! ),..., Qt(sd)) T which yields the same dependence function My as the 
one in (|5|). 



Exhaustiveness of the extremal t dependence structure 

In general, regular variation of the radial variable Rd ensures the MDA condi- 
tion for an elliptical distribution and is a necessary an d sufficient condition for t he 
presence of asymptotic dependence (cf. Theorem 4.3 in lHult and Lindskod (|2002l )). 
Since we can write the (zero-truncated) components of an elliptical random vector 
as Rd(AU)J (j — 1, d), these are also regularly varying in the case of Rd & RV Q 
due to Breiman's theorem. The value of the location vector \i is irrelevant in this 
context. For the other elliptical distributions in the Gumbel or Weibull MDA, as 
for instance the multivariate normal distribution (Gumbel), one always obtains in- 
dependence in the max-stable limit G. If Rd £ RV Q , we can use Proposition 12.11 
and write M(z) as 

M(z) = M a , E . = {E[([/+) a ]} _1 E max {^[(^U)/]"} 

j — 1 , . . . , d 
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with A*(A*) T = £*, and we see that the radial variable enters the dependence 
structure only through its index of regular variation. Since the multivariate t distri- 
bution covers the full range of tail indices (equal to the general degree of freedom) 
and correlation matrices £*, the extremal t dependence structure is exhaustive 
within the class of asymptotically dependent elliptical distributions. 



3. A spectral construction of the extremal-f process 

We first provide a multivariate spectral construction for the extremal t distri- 
bution based on elliptical distributions. 

Theorem 3.1 (Multivariate spectral construction). Suppose the following items 
are given: 

• a tail index a > 0, 

• iid replications Xi of an elliptically distributed random vector X = (X\, X4) 
with dispersion matrix E* (a*^ = 1 for j — 1, d) and location vector fi = 
such that the expectation m a = E[(X^ ) a ] is non-null and finite and 

• a Poisson process {Vi} ~ PIUs/L(av~( a+1 ^ dv) on (0, oo). 
Define the componentwise maximum 

Z = m-^ 1 max T^X; . (7) 

i=l,2,... 

Then Z follows the extremal t distribution with a-Frechet marginal distributions 
and dependence function M a ^* . 

Proof. Due to infinite number of points in the Poisson process {Vi}, we can re- 
place X by X + in the construction 0. By taking Z to the power of a, i.e. 
m^ 1 maxj=i ) 2,... VtO^t ) a ■> we obtain a special case of the construction (j4j which 
proves the max-stability and the a-Frechet marginal distributions of Z. Since 
Vi ~ $ a , the radial variable in the elliptical random vector ViXi is regularly 
varying with index a due to the variant of Breiman's theorem from Lemma 2.3 in 



Davis and Mikoschl (|2008l ). Thus ViXi is in the MDA of the extremal t distribution 
and at the same time in the MDA of Z according to Corollary 12.11 We conclude 
that Z follows the extremal t distribution with dependence function ([6]). 

As a direct application of Theorem l3.11 we are now able to present one possible spec- 
tral representation of extremal t processes via the corresponding Gaussian process. 



Theorem 3.2 (Spectral representation of extremal t processes). Suppose the 
following items are given: 

• a tail index a > 0, 

• iid replications W; of a standard Gaussian random field W on S C W with 
correlation function Cov* and 

• a Poisson process {V,} - PRM(av-( a+1) dv) on (0,oo). 
Then the process defined by 

Z = {Z(s)} = |m Q - Q_1 max Wi(s) j (8) 

with m a = ^Jlx 1 2°' 5 ("~ 2 )r(0.5(a: + 1)) and T(-) the Gamma function, is an ex- 
tremal t process with a-Frechet marginal distributions. Its dependence structure is 
characterized by a general degrees of freedom and the dispersion function Cov* . 



G 



Proof. It remains to verify the value of m a . Using the variable transformation 
y = 0.5a: 2 yields 

' />oo 

///„ = / a; Q (27r)- a5 cxp(-0.5a; 2 )(ia; = / (2y) a5Q (27r)-°- 5 exp{-y)(2y)~ - 5 dy , 
la J o 

and gathering the involved constants leads to the desired representation of m a . 



4. Discussion 



Taking {Z(s)} in ([8]) to the power Z a (s) establishes unit Frechet marginal distri- 
butions as in the construction ( ffl with Q = W" . Clearly, we identify the extremal 
Gaussian process for a = 1 (cf. lSchlatherl (|2002h 'l. 



Remark 4.1 (Extremal coefficient). When d = 2 and a\ 2 — in Theorem 
\3.1l the range of the extremal coefficient covers the open interval (1.5,2): As the 
degree of freedom v — a tends to infinity in the univariate t- distribution 

converges towards the normal distribution and its variance tends to such that 
M(l) = 21im„^ 00 i I , + i(l | 0, (1 + i/)" 1 ) = 2 in ©. As i> tends to 0, we observe 
M(l) = 2 lim^oo ^+1(1 | 0, (l + v)- 1 ) = 2ii(l | 0,1) = 2( 7 r- 1 arctan(l) +0.5) = 
1.5. This helps understand the long-range dependence structure in models for ex- 
tremal t processes since the applied correlation functions are usually positive and 
approach as the distance between two points increases to infinity. 

In particular, extremal t processes can be considered more flexible than extremal 
Gaussian proc e sses o r Brown-Resnick processes which are both special cases; see 
Davison et alJ (|2012h for the case of the Brown-Resnick process which arises for 
some correlation structures as a tends to infinity. Moreover, the formulation of the 
dependence function |6]) for the extremal G aussian process (y = 1) is more general 
than the bivariate expressions obtained by ISchlather ( 20021 ) and lends itself more 
easily to interpr etation. Combining extrema l t mo dels with a random set approach 
as described in Davison and Gholamrezaeel ( 2012 ) further enlargens the range of 
available models and dependence structures. 

Remark 4.2 (Simulation). Theorems \3.1\ and \3.2\ allow to simulate extremal t 



distrib utions and processes with the method devised in Theorem 4 of \Schlathe7 



(2002). When the degree of freedom is large, the computational complexity can be- 
come very restrictive, making it diffi cult to assure a go o d qua lit y of simula t ion. H ow- 



r, the Hiisler-Reiss distribution !(Husler and Reis'i (198& ); Falk et al\ \20l\ )) 



in 



the multivariate case and the Brown-Resnick process in the infinite- di mensional case 
are ad equate proxies for some dispersion structures, as explained in Davison et all 

(MR). 



Future research should explore in more detail up to which general degree of freedom 
a the simulation procedure is numerically feasible, and if the Brown-Resnick process 
provides an adequate substitute around and beyond the "critical" value of a. The 
spectral construction ([8]) further opens t he way for tackling condi t ional simulation in 
th e theoretical framewo rk de veloped b y Dombrv and Eyi- Minkol(|201ll ) and applied 
in lDombrv et all (|201ll ) and lDombrv and Ribated (|2012t ). 
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